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I. INTRODUCTION 

In the last few years the evidence that the Universe 
contains a large proportion of some not ordinary stuff 
which makes it expand acceleratedly is getting stronger 
grounds (see Q for a recent review). Due to such con- 
siderations, quite a few names have been added to the 
original list of contents of the Universe along its differ- 
ent epochs. Among those, quintessence and phantom 
energy are the most popular, and, particularly the lat- 
ter, has triggered a feverish activity in many directions, 
one of them being the investigation of the unusual geo- 
metric properties of the cosmological models they lead to. 
Phantom universes have rather awkward singularities 
which go into the basket of "cosmological milestones" 
along with other geometric curiosities such as bounces or 
extremality events appearing in other sorts of universes. 

Observationally, phantom universes seem to be pre- 
ferred over geometry- wise dull LCDM cosmologies 
thus making legitimate the sort of questions addressed 
here about the fate of phantom universes. We are go- 
ing to perform an innovative analysis of those models 
in conjunction with all the other FLRW models in the 
literature, which will bring some surprises to build on 
the atypicality of phantom cosmologies. In addition, our 
analysis opens new paths for the exploration of other sort 
of milestones, such as sudden singularities which have re- 
ceived considerable attention recently. 

We have borrowed the denomination "cosmological 
milestones" from Q, where FLRW cosmological models 
were analysed in terms of a generalised power expansion 
of the scale factor. The appearance of polynomial scalar 
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curvature singularities and derivative curvature singular- 
ities, together with the satisfaction of energy conditions 
were shown to depend most generally on just the first 
three terms of the expansion. 

Clearly, there are many interesting geometrical fea- 
tures which are elusive to studies of that sort. Since 
the usual definitions of singularities are related not only 
to the properties of the curvature tensor but also to the 
existence of causal geodesies that cannot be extended to 
arbitrary values of their proper time (geodesic incom- 
pleteness) [j| or even of general causal curves with the 
same property (b- incompleteness) 0, Q, it is of inter- 
est to analyse singularities in general FLRW cosmologies 
within this framework, as it was done in for sudden 
singularities @. It is relevant to do so because causal 
geodesies describe the trajectories of observers subject 
just to gravitational forces. Note that curvature is a 
static concept, in the sense that it only reflects what 
happens at each event, whereas features derived from 
tracking the observer along its trajectory are more dy- 
namical, and somewhat more enlightening. Thus, our 
study covers key issues that were overlooked in recent 
related classifications 0, El 

We begin therefore in section [n] by arranging geodesic 
equations for FLRW cosmological models in a suitable 
fashion for integration in terms of a generalised power ex- 
pansion in coordinate time. We proceed then in section 
IIIII to analyse the behavior of lightlike geodesies in these 
models, which sets the foundations allowing to check 
whether the singularities that are found are strong or not 
according to the usual definitions reviewed and refined in 
section lTVI Timelike geodesies are integrated in sectionfVl 
and the strength of their singularities is dealt with in sec- 
tion ^T] The paper ends with a discussion of the results 
in section IVHI Special remarks are done throughout the 
paper regarding observationally allowed/favoured phan- 
tom cosmologies, because there is a peculiar class of such 



2 



cosmologies which persist to stand out of the crowd of all 
phantom models, as long as their geometrical properties 
are concerned. 



II. GEODESICS IN FLRW COSMOLOGICAL 
MODELS 



We consider spacetimes endowed with a FLRW metric 
of the form 

ds 2 = -dt 2 + a 2 (t) {f{r)dr 2 + r 2 (dO 2 + sin 2 9d(f> 2 )} 



f(r) = 



1 — kr 



;, k = Q,±l, 



(1) 



As in we assume that at a coordinate time to, a 
singular event or cosmological milestone comes up in such 
spacetime. To allow our results to be most general, we 
just require the scale factor a(t) to have a generalised 
Puiseux expansion around the event at to, 

a(t) = c \t - t r° + Cl \t - + ■ ■ ■ , (2) 

where the exponents rji are real and ordered, 

Vo <Vi < ■■■ 

This framework covers every proposal of FLRW cos- 
mological model in the literature. 

At first order, a model admitting such expansion be- 
haves as a power- law model of exponent 770, which in the 
case of a flat universe, k = 0, would correspond to a 
linear equation of state for the cosmological fluid, 

p = wp , w = — 1 + 2/3?7o • 

In order to have a positive expansion factor, we require 
Co to be positive. 

Depending on the value of 770, several types of cosmo- 
logical milestones arise Q: 

• t]q > 0: the scale factor vanishes at to and generi- 
cally we have a Big Bang or Big Crunch singularity. 

• 770 = 0: the scale factor is finite at to. If a(t) is 
analytical, the event at to is regular. Otherwise a 
weak or sudden singularity appears [fj . 

• 170 < 0: the scale factor diverges at to and a Big 
Rip singularity appears. 

Since the singular event at to is approached from just 
one side (the past for a Big Bang, the future for a Big 
Crunch singularity), there is usually no need to consider 
absolute values in the expansion (J2J, except in the case, 
for instance, of sudden or weak singularities 9], which 
have been seen to be traversable [8| for geodesic ob- 
servers. 

In order to avoid dealing with signs, we consider singu- 
larities in the past, t > to- Of course the analysis is valid 
also for singularities in the future, since the equations are 



time-symmetric, and occasionally we will comment what 
would happen if the singularity lies in the past. 

We consider causal geodesies parametrised by their 
proper time r, (t(r), r(r), 9(t), 4>(t)) . This means that 
the velocity u of the parametrization (t, r, 9, cf>) satisfies 

S = —gijx'x- 1 , x 1 ,x J — t,r,9, <\> , (3) 

where S takes the zero value for lightlike geodesies and 
the value one for timelike geodesies. It takes the value 
minus one for spacelike geodesies, but since we need just 
causal curves for our analysis, we will discard the d = — 1 
case. The dot stands for derivation with respect to proper 
time. 

Condition defines proper time up to a change of 
scale and a translation, f = At + B, and therefore the 
parametrization is also called affine parametrization. 

Geodesic equations are quasilinear in the accelera- 
tion (t, r, 9, (f>) and depend on the metric components 
through the Christoffel symbols, 

(4) 



T i jk x j x k = 



1 



T )k = 7}9 d {9lj,k + 9lk,j - 9jk,l} ■ 
For a FLRW cosmology they may be written as 



1 - kr 2 



(5) 



r z - aa!r z 9 Z + sin^ 9<j> 2 (6a) 



_ a 1 ■ . kr . 2 

r = -2—tr - — -r 

a 1 — kr z 



+ (1 - kr 2 )r 



-2—i9 - -i-9 + sin 9 cos 6<j) 2 
a r 

a' ■ ■ 2 • 

-2— t(j) fc/) - 2 cot 

a r 



(6b) 
(6c) 

(6d) 



where the comma stands for derivation with respect to 
coordinate time t. 

Taking into account that orbits of geodesies in spher- 
ically symmetric spacetimes remain in equatorial hyper- 
surfaces, they can be fit in a hypersurface 9 = tt/2 by 
choosing coordinates accordingly. This allows a simpli- 
cation of the system of equations, 



t 



1 - k 



rf — aa r 



'^2:2 



(7a) 



kr 



2 i' 

a 1 — kr 2 

2 -; 

-2— td> f(p . 

a r 



r 2 + (1 - kr 2 )rct) 2 (7b) 
(7c) 



Finally, due to the existence of isometries, the following 
conserved quantities of geodesic motion exist: 



Pi = a(t)^ f(r) cos 4>r- jj-^ sin cfxf) } (8a) 



P2 = a{t)\ f(r) sin + — cos cjxj)} (8b) 



L = a{t)r 2 



(8c) 
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They allow us to reduce (J to a simple set of first order 
differential equations: 



e = s- 



P 2 + kL 2 
a 2 (t) 



Pi cos 4> + Pi sin < 



a 2 (t)f(r) 



L 



a 2 {t)i 



(9a) 
(9b) 
(9c) 



These constants of geodesic motion are related to an- 
gular momentum and linear momentum, 



p2 = p2 + p2 ^ 



(10) 



and even allow us to obtain the equation of the orbits, 
namely, 



f(r)r 



L 



P2 cos (f) — Pi sin <j> 



(11) 



which is just the equation of a straight line in polar co- 
ordinates in the flat case k = 0: 

Pix ~P lV = L, 

as it was to be expected. 

It may be seen that the distinction between differ- 
ent types of universes (flat, open, closed) appears only 
through the constant k in (|9a|) and in the function f(r) 
in (|9b|) . This function may be factored away in that 
equation just by taking 



R 



Pi cos + P 2 sin ( 



arcsinhr k = — 1 
R = { r k = 

arcsinr k = 1 



Therefore, the relevant information for geodesies is en- 
coded in the scale factor a(t). We make use of this fact 
on analysing the behaviour of causal geodesies. 

From © we learn that just the equation for t needs to 
be solved, since the equations for the other coordinates 
are reduced to quadratures once the solution of (|9a|) is 
known. 

We can forget the equation for </> along a geodesic since, 
due to homogeneity and isotropy of the FLRW universe, 
the origin may be located at any point on the geodesic 
and hence this appears as a straight line with zero angular 
velocity: 



(12a) 
(12b) 




a 2 (t)f(r) • 



Possibly a quicker way to reach this result is consid- 
ering beforehand that geodesies are straight lines due to 
the homogeneity and isotropy of the spacetime and that 



9r = d r / f(r) is a generator of an isometry along one of 
these lines. Hence, 



±P = u- 



W) 



a 2 (t)f(r)f 



is a conserved quantity of geodesic motion. The equation 
for t is derived from the normalization condition @, 

6 = i 2 -a 2 (t)f 2 (r)r 2 . 

We consider future-pointing geodesies and therefore we 
take i > 0. 

Now we may begin to draw information about causal 
geodesies from their equations. Following p| we take 
causal geodesic completeness as a minimum condition for 
a spacetime to be considered singularity-free. Therefore, 
we analyse the cases where causal geodesies are incom- 
plete, that is, where they cannot be extended to arbitrar- 
ily large values of their proper time r. 

However, we must bear in mind that, since we have 
no guarantee that the coordinate chart that allows us to 
write the metric in the form covers the whole uni- 
verse, some conclusions about incompleteness may not 
be correct if the spacetime is extendible to a larger one. 
That is, a geodesic may leave the portion of spacetime 
depicted by our coordinates in finite proper time, but not 
the universe itself. Therefore, some of the singularities 
we may encounter may not be real, since the universe can 
be extended. 

This is the case of Milne universe, which is the case 
of for k = — 1, a(t) = t. This universe can be re- 
duced to a portion of empty Minkowski spacetime by the 
coordinate transformation 



T = t\fl 



R = rt 



which covers just the region inside the null cone T 2 = R 2 . 
It is, therefore, a geodesically complete and singularity- 
free spacetime, but it appears singular in the Milne form, 
since it can be extended to the whole Minkowski space- 
time. 



III. LIGHTLIKE GEODESICS 

The lightlike case is fairly simple and can be explicitly 
integrated for the time coordinate, 



a{t)i = P 



a(t') d£ = P{t - t ) 



Close to to, the leading term in the power expansion 
of equation l|12a|) is the one with lowest exponent, 770- I n 
many cases, in order to analyse the singular behavior of 
geodesies near io ; w e just require the first term of the 
power expansion, 



a(t) ~ co|* — *o| 



</(> 
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which provides the time coordinate at first order in terms 
of proper time r, after integrating (|12a|l . 

f c |t' - ioT dt' ~ P(r - T ) , 

•/tn 



CO 



(r - r ) 1/(1+ " o) , (13) 



for 770 =/= — 1. If 770 = — 1, the leading term is exponential, 
t ~ < + Ce Pr/co . 

Other cases which require a different treatment, involv- 
ing more terms of the power expansion, are those with 
770 = (1 — n)/n with n a positive natural number. We 
elaborate further on those cases below. 

From either the expression of t in 1|13|) or t and its 
derivatives one gets 



t 
i 



P 

co 



|i-f |-"° 



CO 

VoP 2 



\t-to\ 



\t - hr 2 ^- 1 



A„|t-t | 1 - n -"' ?0 



(14a) 

(14b) 
(14c) 



with 



A„ = (-l)" +i 77o 



((n - 1)770 + n - 2) f^j . (15) 



Several possibilities arise, since (|14c|) implies that if 



Wo < 



1 — n 



hn — 3 
3(1 - n) 



<w<-l 



(16) 



there is no blow up in any of the derivatives of order 
lower or equal than n. The latter condition is not too 
stringent, that is, such cases appear often, so it is of 
utmost relevance to get further insight into the details of 
the different subcases one can distinguish: 

• 7/0 > 0: This case includes all classical matter con- 
tents (for flat universes, scalar field 770 = 1/3, radi- 
ation 770 = 1/2, dust 770 = 2/3. . . , with w > —1). 
Since the exponent 1/(1 + 770) is lower than one, t 
is not differentiable at to and the derivative t blows 
up. 

• 7/0 G (—1/2, 0): It corresponds to w < — 7/3 for flat 
power-law models. In this case i does not blow up 
at to, but t does. 

• 770 G (—2/3,-1/2): It corresponds to w G 
(—7/3, —2) for flat power-law models. In this case 
i does not blow up at to, but i does. 



•Wo € (^T i '^r) : I* corresponds to w G 
( 3{2-n) ' 3^i-k) ) for nat P ower -law models. The 
derivative t" -1 ) does not blow up at to but t n ^ does. 

• 770 < —1: It corresponds to w G (—5/3,-1) for 
power-law models. According to (JT3J, the time co- 
ordinate along the lightlike geodesic is dominated 
by a negative power of proper time r. This means 
that these geodesies never reach to, since t — to only 
vanishes when r tends to infinity. Typically, best 
fit phantom models have a value of w within this 
range, so it seem very likely that if the universe is 
phantom its geodesies are going to have this pecu- 
liar behavior. 

The limit cases where 770 is of the form (1 — n)/n (0,- 
1/2,-2/3, -3/4, ...,-1) fall out of this classification since 
the derivative t n+1 ^ vanishes, as it follows from l|15|) . If 
every exponent 77^ in the generalised power expansion of 
a(t) were of this form, none of the derivatives of the 
time coordinate along lightlike geodesies would blow up. 
These are extremely fine-tuned cases, so we will not con- 
sider them any further, and we will then turn back to the 
cases 77, = (1 — n)/n for i = only. 

For this analysis we have to resort to the next term of 
the power expansion, c\(t — to) Vl , but if the term does 
not provide sufficient information, one would have to keep 
adding terms till a satisfactory expression is obtained. 

Up to second order, 



t = 



P P 



PC A 



a(t) co 



-\t-to\- r > -~-^\t-to\ r >i- 2r >° 



r 2 



we see that, after the contribution of the first term to the 
derivative t n+1 ^ vanishes, 

- |i-t | I ' 1 - ( ™+ 2 ' , "'-' i , 

the leading term for lightlike geodesic behavior in these 
cases is the one with 771. 

Let us take a look to the new cases arising: 

• 770 = 0. In this case the scale factor is neither zero 
nor it tends to infinity at to and, in principle, equa- 
tions Ijl2a|) are regular at to as it was pointed out in 
|j| . This is the case of the models proposed by Bar- 
row . How regular/singular geodesies are in these 
spacetimes depends on the next exponent, 771 > 0, 
in the expansion of the scale factor. For any such 
value t does not blow up at to, but some remarks 
are in order. 

If 77! G (n — 2, n — l), n > 1, the derivative t n_1 ) 
does not blow up at to but t n ' does. The phe- 
nomenon called sudden sing ularity 0, which has 
aroused much interest (3, [HI U2, 1 > corresponds 



In |IJ and some of the references inspired by that work llll Il4l , 
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to the models with rji > 1, which have therefore 
nonsingular t at to, and were shown to have well- 
defined geodesies around to in || ■ 

If rji is a natural number, the same reasoning is ap- 
plied to the first exponent i]n which is not natural. 
If every exponent r\i is natural, then no derivative 
of lightlike geodesies diverges at to- This is the 
case of course of a non- vanishing analytical a(t) 
in the vicinity of to, such as, for instance, in de 
Sitter, a(t) = cosht, and anti-de Sitter universes, 
a(t) = cost. 

• rjo = —1/2. For all of these cases t does not blow 
up at to- Since f ~ |t — to\ m , the third derivative 
blows up if rji is negative. 

Since t™> ~ \t - t |' 7l+(3_n)/2 in this case, if rfr G 
({n - 4)/2, (n - 3)/2), n > 2, the derivative t™" 1 ) 
does not blow up at to but t™- 1 does. 



Vo 


m 


i 


i 


"t 




(0,oo) 


(Vo,oo) 


oo 


oo 


oo 


OO 





(0,1) 


finite 


oo 


oo 


oo 


(1,2) 


finite 


finite 


oo 


oo 


(2,3) 


finite 


finite 


finite 


oo 


(-1/2,0) 


(»?o,oo) 


finite 


oo 


oo 


oo 


-1/2 


(-1/2,0) 


finite 


finite 


oo 


oo 


(0, 1/2) 


finite 


finite 


finite 


oo 


(-2/3,-1/2) 


(??o,oo) 


finite 


finite 


oo 


oo 


-2/3 


(-2/3,-1/3) 


finite 


finite 


finite 


oo 


(1-n 2-n\ 
I n ' n— 1 ) 


(Vo,co) 


finite 


finite 


finite 


oo 


(-oo,-l] 




/ 


/ 


/ 


/ 



TABLE I: Derivatives of lightlike geodesies at to. A slash 
indicates the cases where to is never reached. 



IV. STRENGTH OF SINGULARITIES ALONG 
LIGHTLIKE GEODESICS 



• tjo = —2/3. In these models t is finite at to- Since 
i ~ \t — t | ,,1+1 / 3 , the fourth derivative blows up 
if 7/1 < -1/3. 

Since t") - |t-t |' n+(:5 ~ n)/3 , if r?i G ((n-6)/3,(n- 
5)/3), n > 3, the derivative t"^ 1 ) does not blow up 
at t but t n > does. 

Therefore we see different levels of singular behavior 
at t = to along lightlike geodesies. As tjo decreases to 
get closer to the value —1, the regularity of the geodesic 
improves in the sense that one has to go to derivatives of 
higher order to find singular behavior. This limit model 
has been named superphantom [l5j | in the case of flat 
power- law cosmologies with w = —5/3. For 770 greater 
or equal than —1 these geodesies do not even reach the 
event at t . 

As we have already pointed out, the geodesic equation 
for r does not add any further information, as we see in 
the equation of the orbit 



dR 
~dt 



R 



t 



_L_ ^ |t-to|-"° 
a(t) Co 



which is integrable close to to if 970 < 1- That is, it adds 
no new information, since cases with 770 > 1 have already 
been shown to be singular. 



the treatment is purely phenomenological and the focus is on a 
quest for ad hoc para metr izations of a(t) leading to such singular 
behavior to arise (see Il3l for an approximate reconstruction of 
the equation of state or Il4l for statefinder parameters). Inter- 
estingly, sudden singularities can appear in more solid contexts, 
for instance in some braneworld models 



As we have shown in the previous section, the 
"strength" of singularities decreases qualitatively as the 
exponent 770 decreases. It would be interesting to check 
these differences in behavior with the usual definitions of 
strong singularities. 

The idea of a strong singularity was first introduced 
by Ellis and Schmidt |17j . A singularity is meant to be 
strong if tidal forces exert a severe disruption on finite ob- 
jects falling into it. There have been several attempts to 
provide a rigorous mathematical definition for this idea. 

The finite volume is considered to be spanned by three 
orthogonal Jacobi fields which form an orthonormal basis 
with the velocity of the incomplete geodesic. According 
to Tipler |18| , the singularity is strong if the volume tends 
to zero as the geodesic approaches the value of prop er 
time where it meets its end. Krolak's definition |l9| is 
less restrictive, since it just requires the derivative of the 
volume with respect to proper time to be negative. This 
definition has been further refined in 20] . 

However, these definitions are meant for focusing grav- 
itational forces, that is, RijU l u J must be non-negative for 
timelike and lightlike observers with velocity u. Hence, 
these definitions leave out the possibility of Big Rip sin- 
gularities, but they can be extended to these cases just re- 



versing signs. For instance, for negative Ri 



Tipler's 



definition requires a volume tending to infinity as the 
geodesic meets its end and Krolak's definition requires 
positive derivative of the volume close to the end of the 
geodesic. 

Fortunately, one is not to construct the basis of Jacobi 
fields in order to check these definitions. There are nec- 
essary and sufficient conditions due to Clarke and Krolak 
pl| related to integrals of Riemann tensor components 
along the incomplete geodesic. They are not affected 
by the inclusion of Big Rip singularities since they just 
require components of the Riemann tensor to blow up 
along causal geodesies. In order to apply them we need 
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the components of the Riemann tensor, 

T3t 

^rtr 



pt 

n ete 



p 



pt 

r 2 sin 2 i 



• For r/o € (— 1, 1), k ^ or ry > — 1 5 k = 0, using 

ecu, 



RijU l u J 



| t _ ^1-2(^+1) 



pi* _ p 

-fU^ — -ft 



7?^ — — 



2% 



(77o + l) 2 |r-r |2 ' 



pr- Rr jn0 nfl 

^flrS = -"^0 _ _£Vr£ = -"^ 

r 2 r 2 sin 2 / 2 r 2 si 



rrcj) 



R 



8H =a' 2 + k 



omitting the ones that can be obtained by the symmetries 
of the tensor, and Ricci tensor components, 



Rt 



3a" R rr _ Ree _ R4 



1*2 cin^ 1 



aa" + 2(a' 2 + k) 



For FLRW models conditions are simpler, since the 
Weyl tensor vanishes and therefore only conditions re- 
lated to the Ricci tensor are relevant. According to |2lj . 
a lightlike geodesic meets a strong singularity, according 
to Tipler's definition, at proper time To if and only if the 
integral of the Ricci tensor 



dr' / ,/r" /,',,,/•,/•' 

I I JO 



(17) 



diverges as r tends to to. 

For Krolak's definition the condition is less restrictive: 
a lightlike geodesic meets a strong singularity at proper 
time tq if and only if the integral 



dr'Rijtfu 3 



diverges as t tends to t . 

In our case, the velocity of the geodesic is 



(18) 



u l = t 



P 



r • 1 

u = r = ± 



fa 2 



and therefore the component of the Ricci tensor measured 
along the geodesic is 



RijU l u j = 2P 2 



a' 2 + k a" 



At first order in our generalised power expansion, we 
have two cases, depending on whether the curvature term 
is leading or not, 



2P 2 m 



2kP 2 



R ^ UJ ^ c 2\ t _ tQ \2 { , + l) > 4\ t -t \^o 



produces a logarithmic divergence with Tipler's 
definition and an inverse power divergence with 
Krolak's for ^ ijo > —1 and therefore we have 
a strong singularity in these cases. 

• There is a subcase left, 770 = 0, for which the ap- 
proximation at first order leaves 



RM * 2F 2 f4 c ^(m-l) 



4 4\t-t Q \ 2 -^ 

2kP 2 2P r > 1 c 1 rj 1 (Ti 1 - 1) 



2-m 



which provides no divergent integral with Tipler's 
definition, since 771 > 0, but provides one with 
Krolak's one if 771 € (0,1). This generalises the 
result of Q, since there it was shown that sudden 
singularities, a special case with T)q = 0, rji > 1, 
were not in fact singularities. 

For rji = 1 we have to use still another term of the 
expansion, 



Rijtfv? ~ 2P 



, I k + c\ c 2 -q 2 {m - 1) 



4 4\t-to\ 2 -*\ 
2{k + c 2 )p 2 _ 2P^c 2 m{m - 1) 

4 Crf+V-Tol 2 "* ' 

which shows that this subcase does not produce a 
divergent integral since 772 > 1. 

• For 770 = 1, 



i?ytiV ~ 2P 



4 + k, 



2 L 



t-tn 



we see that both integrals are divergent, since the 
exponent is smaller than —2, unless k = — 1, cq = 1, 
values for which the Ricci tensor vanishes at first 
order, because at this order it is Milne model. We 
are to resort then to the next term, 77! > 1, 



RijU l u J 



2P 2 m { m - 3)ci 

l^-toM 1 

2 (m-3)/2p{m-i)/2 Tjl ( T}1 _ 3 ) Cl 
|t-t |( 5 -"i)/ 2 



• First, the case 770 < — 1, which is complete, since 
lightlike geodesies never reach t = to, as we have 
seen. 



which produces no divergent integral for Tipler's 
definition but for Krolak's one the singularity is 
strong if 771 < 3. 
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k 


CO 


Tipler 


Krolak 


(-oo,-l] 


(770,00) 


0,±1 


(0,oo) 


Regular 


Regular 


(-1,0) 


Strong 


Strong 





(0,1) 


Weak 


Strong 


[l,oo) 


Weak 


Weak 


(0,1) 


fa, 00) 


Strong 


Strong 


1 


(1,00) 


0,1 


Strong 


Strong 


(1,00) 


-1 


(0,l)U(l,oo) 


Strong 


Strong 


(1,3) 


1 


Weak 


Strong 


[3,oo) 


Weak 


Weak 


(1,00) 


fa, 00) 


0,±1 


(0,oo) 


Strong 


Strong 



TABLE II: Degree of singularity of null geodesies around to. 



We need another term, 772 > 3, to check the regular- 
ity of the 771 = 3 subcase with Krolak's definition, 



RijU l u 3 



2P 2 m(rn -3)ca 



I* — *o| 



5-»72 



and we find that is similar to the 771 contribution. 
Hence it does not diverge for r; 2 > 3. 



• Finally, in the cases with k ^ and 770 > 1 the 
leading term in the Ricci tensor component is the 
curvature one, 



RijU l u J 



2kP 2 



2kP 2 



a 4 c 4 |i-i n | 4 "° 
2k(l + r? )" 4/(1+I ' o) P 2 - z WU+'?o) 

c 4-4WU+'Jo)| T _ ro |4W(l+'7o) ' 

which also provides divergent integrals since the ex- 
ponent of the denominator is larger than 2. Hence, 
the singularities are strong in these cases either. 



We conclude that for 7/0 > — 1 lightlike geodesies in 
all models meet a strong singularity except for the cases 
7/0 = (and 771 > 1 with Krolak's definition) and 770 — 
1, k = —1, Co = 1 (771 > 3 with Krolak's definition), 
which are regular. These are the only cases, together 
with 770 < —1, that are null geodesically complete, even 
though the curvature is singular also for these models, as 
it was shown in 

We notice that the different behavior of geodesies for 
positive and negative rjo does not quite influence the 
strength of the curvature singularity at to. Generically 
models with a Big Rip have null geodesies with deriva- 
tives that do not blow up at to whereas all derivatives of 
null geodesies in models with a Big Bang or Crunch are 
infinite. 



V. TIMELIKE GEODESICS 

For timelike geodesies the relevant equation is, at first 
order of the power expansion, 



p2 

9" 

a 



P 2 

j l + —(t-t )-^o ; 



(19) 



which can be solved explicitly in terms of hypergeometric 
functions, 



1 

V 



(t-t )F ( T (t-t 



1 



p2 



,-2i)o 



II 



TO , 



where F is the hypergeometric function, but we shall not 
use this expression. 

It is clear, as it happened for lightlike geodesies, that 
for 770 > the geodesic is singular at t = t , since t blows 
up there, unless P is zero, which is a trivial regular case, 



t - t = T - T 



r = r , 



the comoving congruence of fluid wordlincs. 

On the contrary, this derivative is well defined and 
takes the value one if 770 is negative. In this case we 
may approximate t in the vicinity of to, 



Z6 



(20) 



In order to carry out the analysis of the behavior of 
these geodesies, we need expressions for higher deriva- 
tives of coordinate time t, 



PV 



PV 



t = P 



P 2 Vo(t - tp)- 2 ^- 1 

r 2 
c 

2 f3a' 2 a"\ 

3" 

a 1 a 3 



with 



\i — 



_ P 2 77o(277o + l)(f,-£o)- 2 "°- 2 
~ r 2 

' ^ \n(t-t )- 2 »°- n+1 . 

(-l)" +1 P 2 277 ---(277o + n-2) 
2c 2 



(21a) 
(21b) 
(21c) 



(21d) 



From these expressions we may draw valuable informa- 
tion about timclike geodesies around to'- 

• 770 > 0: As it happened in the lightlike case, time- 
like geodesies are singular at to, since the derivative 
t blows up. For flat power-law models it corre- 
sponds to w > — 1. 
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• r?o S (—1/2,0): Again as it happened for lightlike 
geodesies, the derivative i blows up at t whereas 
i does not. It corresponds to w < — 7/3 for flat 
power-law models. 

• r]o G (—1,-1/2): For these cases we find the first 
difference with lightlike geodesies. The derivative 
i blows up at to, but t does not. They correspond 
to flat power-law models with w € (—7/3, —5/3) . 

• Hi 6 (^ir^^ir); n > I: Whereas lightlike 
geodesies did not reach to in finite proper time 
in these models, timelike geodesies do, with reg- 
ular t n ~V but with infinite t n) at t . The corre- 
sponding flat power-law coefficients would be w € 

G 



3n-2 3n+l 
3(2-n) ' 3(l-n) 



Again for the limit cases r/o = =^ (0, —1/2, —1 . . .) the 
contribution of the term with exponent 770 to the deriva- 
tive vanishes and we have to resort to the next 
term in the expansion with non-vanishing contribution 
to higher derivatives, 



f 



S(*-*o)- 4 " 



P 2 Ci 



(t - to)" 1-8 " 



+ 



(22) 



which is the term with exponent 771 — 3770 and therefore 
the relevant contribution to the derivative t n+1 ' is of the 
form (t-t ) m - 3v °- n . 

Let us analyse some of these cases: 

• 770 = < 771 : The discussion is entirely similar to 
the one for lightlike geodesies. The scale factor does 
not vanish at to and therefore these are sudden or 
weak singularities. Since t ~ \t — to| r?1 ) they have 
finite t at to and the derivative t is finite for 7/1 > 1. 
If 7/1 € On — 2, n — l) , n > 2, the derivative t"^ 1 ' is 
also finite at to but t n ^ is not. 

Again, if 771 is natural, we would have to resort to 
the first exponent rjjq which is not natural and if all 
of them are natural, then no derivative of timelike 
geodesies diverges at to. 

• rjo = — 1/2. In these cases t is finite at to. Since 
i ~ \t — to| m-1 ^ 2 , the third derivative is finite if 
771 > 1/2. When 771 e (77 - 7/2,77 - 5/2), n > 3, 
the derivative i™ -1 ) is finite at to and i™- 1 is not. 

• rjo = —1. Now i is finite at to and t ~ \t — to\ m . 
Hence the fourth derivative is finite if 771 is positive. 
For 7716 (n — 5, 77 — 4) , n > 4, the derivative t n ~ x ^ 

is finite at to whereas t n ^ is not. 

Summarizing, geodesic behavior is similar for timelike 
and lightlike geodesies in models with 7/0 > —1/2, but 
there is a different pattern for the rest of the models. 



770 


771 


i 


i 


"t 




(0, 00) 


(»7o,oo) 


00 


00 


00 


OO 





(0,1) 


finite 


00 


00 


OO 


(1,2) 


finite 


finite 


00 


OO 


(2,3) 


finite 


finite 


finite 


OO 


(-1/2,0) 


(770,00) 


finite 


00 


00 


OO 


-1/2 


(-1/2,1/2) 


finite 


finite 


00 


OO 


(1/2,3/2) 


finite 


finite 


finite 


00 


(-1,-1/2) 


(??o,oo) 


finite 


finite 


00 


00 


-1 


(-1,0) 


finite 


finite 


finite 


OO 


(l—n 2-n\ 




finite 


finite 


finite 


OO 



TABLE III: Derivatives of timelike geodesies at tg 



Differentiability of timelike geodesies improves as 770 de- 
creases, but there are only isolated cases for which they 
are completely regular and this makes a difference with 
the lightlike case. There are no timelike geodesies which 
take an infinite proper time to reach to, as it happens 
with null geodesies with 770 < —1. 

The equation (|12b(l for r does not add further infor- 
mation on the behavior of timelike geodesies either. We 
may tackle the equation of the orbit of the geodesies, 



dR 



dt 



R 



P 



i a(t)^P 2 +a 2 (t) ' 



in a similar fashion. 

If 770 > 0, we get, close to to, 



dR 



dt 



< 



1 

W) 



1 

Co 



is not integrable if 770 > 1, but these are all already sin- 
gular cases. 

If ?7o < 0, we get, close to to, 



dR 



dt 



< 



P 



p 

~2 \t — to 

4 



-2770 



which is integrable for 770 < 1/2. 

Therefore, no new singular behavior appears on consid- 
ering the geodesic equation for r. The radial coordinate 
is singular where t is already singular. 



VI. STRENGTH OF SINGULARITIES ALONG 
TIMELIKE GEODESICS 

Again, it would be quite interesting to know whether 
the singularities encountered by timelike geodesies are 
strong or not according to the usual definitions. 

Conditions like I|17ll8fl are not so simple for timelike 
geodesies, since there are no both necessary and sufficient 
conditions in this case. Those conditions become just 
sufficient if the Weyl tensor vanishes. 
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A timelike geodesic meets a strong singularity, accord- 
ing to Tipler's definition, at proper time To if the integral 
of the Ricci tensor 



dr' 
n Jo 



dT"R i jV?vP 



(23) 



diverges as r tends to t . 

In contrast, for Krolak's definition, a timelike geodesic 
meets a strong singularity at proper time Tq if the integral 



dr' RijU l u J 



(24) 



diverges as r tends to r . 

First we use the comoving fluid worldline congruence, 
with velocity u — dt- In this case, t = 1, proper time and 
coordinate time are the same t — to = r — tq. The com- 
ponent of the Ricci tensor measured by observers along 
this congruence, 



RijU l u 3 = — 



3a" 
a 



3?7o(77o - 1) _ 3770(770 - 1) 
\t-t \ 2 " |r-ro| 2 



produces a logarithmic divergence with Tipler's defini- 
tion and an inverse power divergence with Krolak's one, 
so we may conclude that the singularities are strong for 
all models with 1 =^ 770 7^ 0. Therefore we have: 

• For 1 ^ )/o ^ the geodesies in the fluid congruence 
meet a strong singularity at to. 

• For 770 = 0, we need another term in the expansion, 

3a" ^ 3ci77i(77i - 1) _ 3ci77i(?7i - 1) 
a ~ c \t - tol 2 -" 1 ~ c \t-t \ 2 -^ ' 

and we see that in these cases the integrals do not 
diverge with Tipler's definition, but they do with 
Krolak's for r}\ <E (0, 1). 

• For 770 = 0, rji = 1, we still need another term, 

3c 2 ?72 (772 - 1) 3c 2 77 2 (?72 - 1) 



3a" 
a 



c |t - tol 2 "" 2 c |t-t | 2 -" 2 



in order to check that integrals do not diverge for 
these models with both definitions, since 772 > 1. 

• For 770 = 1, we resort to the second term in the 
expansion, 

3a" ^ 3ci77i (771 - 1) _ 3ci77i(t7i - 1) 
a ~ col* - tol 3 "" 1 ~ c |t-t | 3 -"i ' 

and we find that for these models the integrals do 
not diverge with Tipler's definition, but they do 
with Krolak's for 771 € (1,2]. 



770 


771 


Tipler 


Krolak 


(-00,0) 


(770,00) 


Strong 


Strong 





(0,1) 


Weak 


Strong 


[1,00) 


Weak 


Complete 


(0,1) 


(770,00) 


Strong 


Strong 


1 


(1,2] 


Weak 


Strong 


(2,oo) 


Weak 


Weak 


(1,00) 


(770,00) 


Strong 


Strong 



TABLE IV: Degree of singularity of the fluid congruence of 
timelike geodesies around to. 



This result may be further refined using timelike radial 
geodesies, for which 



u l = t 



p2 

a 2 



r = ±- 



P 



3a" 
a 



2P Z 



Taking a look at the geodesic equation (|19|l for t, we 
notice three different possibilities: 

• 770 < 0: Since t ~ 1, t — to — t — tq close to to. 



• 770 = 0: Now i ~ y/l + P 2 /cq = a close to to and 
so t — to — a(r — tq). 

• 770 > 0: For these cases t ~ P/a close to to as for 
lightlike geodesies. 

Accordingly, there are several cases: 

• 770 < 0: At lowest order, the P— dependent terms 
■ l2 ^k 



a ' 



cf 



a 



a" 



ro |207o+l) 



produce no divergent integral with Tipler's def- 
inition but they do with Krolak's one for 770 € 
[—1/2,0), but it does not matter, since the first 
term was already seen to be divergent, as it is the 
same as for the fluid congruence in all these cases. 

• 7;o = 0: The P— dependent term is essentially the 
same as for lightlike geodesies and we reach there- 
fore the same conclusion: these models produce no 
divergent integral with Tipler's definition, but with 
Krolak's one they do if 771 6 (0, 1). The same hap- 
pens with the first term, which is the same as for 
the fluid congruence. 

• 770 = 1: The P— dependent term is the same as 
for lightlike geodesies. Hence, these cases are all 
singular but except maybe for k = — 1, Co = 1. 
Models with k = — 1, cq = 1 and 771 < 3 are singular 
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rjo 


m 


k 


CO 


Tipler 


Krolak 


(-oo,0) 


(Vo, oo) 






Strong 


Strong 





(0,1) 






Weak 


Strong 




[l,oo) 


0,±1 


(0, 00) 


Weak 


Weak 


(0,1) 


(770,00) 






Strong 


Strong 


1 


(1,00) 


0,1 




Strong 


Strong 




(1,00) 




(0,l)U(l,oo) 


Strong 


Strong 




(1,3) 


-1 


1 


Weak 


Strong 




[3,oo) 






Weak 


Weak 


(l,oo) 


(770,00) 


0,±1 


(0,oo) 


Strong 


Strong 



TABLE V: Degree of singularity of radial timelike geodesies 
around to 



with Krolak's definition. On the other hand, the 
first term, 



3a" _ 301771(771 - 1) _ (2P\ * 3c m ( m -l) 



a c \t-t \ 3 -^ ~ \ <■„ 1 
does not diverge. 



C0\T ~ T I 



• 1 7^ ?/o > 0: The P— dependent term for these 
geodesies is the same as for lightlike geodesies and 
therefore it is divergent in all cases, though the first 
term, 



3a" ^ 3r/ (r/ - 1) 
a ~ \t-t \ 2 



c \ 1+ ™ 3r/o(vo - 1) 



^(1 + Vo) 



2 , 
|t - r 1+ *>° 



does not diverge with Tipler's definition and only 
with Krolak's for 770 G (0, 1). 

Therefore we have so far exactly the same models with 
strong singularities as we found for lightlike geodesies 
plus the 770 < — 1 models, which are null, but not timelike, 
geodesically complete. That is, we know that all models 
with 7^ 770 7^ 1 have strong singularities at to- We don't 
know what happens with models with 770 = 0, though 
those with 771 G (0, 1) have strong singularities according 
to Krolak. And the same happens with models with 770 = 
1, k = —1, Co = 1, though those with 771 € (1,3) have 
also strong singularities according to Krolak. 

Since the condition on integrals of the Ricci tensor 
is not also a necessary condition for the appearance of 
strong singularities, we have to check other ways to get 
information about the 770 = and 770 = 1 models. 

For Tipler's definition |2l|, if a causal geodesic with 
velocity u meets a strong singularity, then the integral 



dr 1 
i] Jo 



dr" |i^,« fc tt'| 



(25) 



diverges as r tends to To for some i, j. The compo- 
nents are referred to a parallely transported orthonormal 
frame. 



Again Krolak's definition is less restrictive and just 
requires that the integral 



l kjl 1 



,k„.l\ 



(26) 



diverges as r tends to To for some i, j. 

We begin again with the fluid worldline congruence, 
u l = 1, for which the only non- vanishing components of 
the Riemann tensor, 



m 



a 
a 



produce a necessary condition which is the same as the 
already studied sufficient condition. 

Therefore, geodesies in the fluid worldline congruence 
meet a strong singularity if and only if 1 7^ 770 7^ (or 
770 = 0, 7fi 6 (0, 1) and 770 = 1, 771 € (1,2) with Krolak's 
definition). 

Radial timelike geodesies show more strong singulari- 
ties, as we may see. We complete the orthonormal basis 
formed by u, 



u = tdt + fd r 
and a vector v, 

i 



pi p 
a z a z f 



P 



v = afrdt H -d r = ±—d t 

aj a 



1 



P 2 



-d r , 



adding the corresponding unitary vectors parallel to de 
and d^. The parallel transport requirement is trivially 
satisfied. 

The 9 and 4> components of the Riemann tensor, 



R%iu k u l 



T>u 4-2. 
JX t f) t i 



H6t L 
,11 



R r g r r 



= -—-PU^ 



a' 2 + k 



Ket) tr 



have similar terms as RijU l u° and therefore produce the 
same results as the corresponding sufficient condition. 
Finally, the v components, 



R l kjiv l u k v J u l 



Ri 



P 4 



1 



la 6 / 2 a 2 f 2 

P 2 ( P 2 
1 + 



p2 

a 2 



V/ 2 



a 
a 



provide a term that has already been discussed and 
therefore we may conclude also that sufficient conditions 
for the appearance of strong singularities along timelike 
geodesies are also necessary, as it happened for lightlike 
ones. 
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VII. DISCUSSION 

We have obtained a thorough classification of singular 
events in FLRW cosmological models in terms of the ex- 
ponents of a generalised power expansion of the scale 
factor in coordinate time around a cosmological mile- 
stone at to. The behavior of causal geodesies has been 
obtained in the vicinity of the event. The first differ- 
ence that has been found is that whereas the velocity of 
causal geodesies blows up at Big Bang and Big Crunch 
singularities, it is finite at Big Rip singularities, as well as 
acceleration and other derivatives, depending on the first 
exponent in the expansion, 770- For sudden singularities 
the velocity is finite and the acceleration may be finite or 
not, depending on the next exponent, 771. 

However this difference of regularity between Big 
Bang/Crunch and Big Rip singularities does not prevent 
the strong character of both types of cosmological mile- 
stones with both Tipler and Krolak's definitions of strong 
singularities. There is only a curious feature in Big Rip 
singularities in models with 770 < — 1 (which are precisely 
those favoured by observations): lightlike geodesies do 
not reach the curvature singularity at to in finite proper 
time and therefore these spacetimes are null geodesically 
complete close to the singular event. Hence photons 
never experience Big Rip singularities and the universe 
would last eternally for them. This feature, however, is 
lost on dealing with timelike geodesies, which reach to in 
finite proper time and meet a strong singularity. 

The only models which allow regular behavior close to 
to are those with r] = and with 770 = 1, k = —1, c = 1. 
The latter ones are Milne universes at first order, which 
are essentially Minkowski spacetime after extending the 
model beyond to. The former ones include models with 
a non-vanishing analytical scale factor, such as de Sitter 
universes, and models with sudden singularities, which 
have finite velocity, but non-finite acceleration or higher 
derivatives of the parametrization of the geodesies de- 
pending on the exponent 771. The larger this exponent 
is, the better the properties of the model are. These cos- 
mologies prevent the formation of strong singularities ac- 
cording to Tipler 's definition, which requires the crushing 
to zero or disrupting to infinity of finite volume objects 
evolving along causal geodesies. With Krolak's defini- 
tion, which requires just a positive derivative of volume 
for Big Bang and Big Rip singularities and a negative 



derivative for Big Rip singularities, strong singularities 
are avoided in models with 771 > 1. This definition seems 
more appropriate, since causal geodesies in models with 
7/1 € (0, 1) do not have finite acceleration and therefore 
geodesic equations would be singular at to, though the 
curves may be extend beyond that event. 

We may compare these results with those studied by 
Cattoen and Visser in where just singularities in cur- 
vature were considered, without taking into account their 
strength nor the behavior of causal geodesies. Those au- 
thors found that the only models without polynomial 
curvature singularities are those with 770 = 0, 771 > 2 
or 771 = 1, 772 > 2, and those with 770 = 1, k = — 1, cq = 1, 
771 > 3. Dealing with derivative curvature singularities, 
the list reduces to models with 770 = and natural expo- 
nents T]i, i > 1 and those with 770 = 1, = — 1, cq = 1 
and natural exponents r\i > 3, % > 1. Derivative cur- 
vature singularities are not reflected in our classification 
since derivatives of the Riemann tensor appear neither 
in geodesic equations nor in Jacobi equations. The ap- 
parent discrepancy between our results and the presence 
of polynomial curvature singularities lies on the fact that 
either geodesies do not reach that singularity or that they 
reach it, but the curvature growth is not enough to form 
a strong singularity. 

Finally, another consequence is that singularities ap- 
pear just in models with vanishing, divergent or non- 
smooth scale factors. From the mathematical point of 
view at least it is worth mentioning that regular models 
are an open set within the family of smooth homogeneous 
and isotropic spacetimes, as it happened for instance with 
inhomogeneous scalar field Abelian diagonal G2 models 
|22| . On the contrary, singular models are not an open 
set, since the vanishing requirement is not generic. 
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